Microscopic calculation of the wobbling excitations by using the Woods-Saxon 

potential as a nuclear mean-field 
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The wobbling excitations of the triaxial superdeformed (TSD) bands in the Lu and Hf region 
are studied by the microscopic framework of the cranked mean-field and the random-phase ap- 
proximation (RPA). In contrast to the previous works, where the Nilsson potential was used, the 
more realistic Woods-Saxon potential is employed as a nuclear mean-field. The wobbling-like RPA 
solutions have been found systematically in the nuclei studied and their characteristic properties 
are investigated in details. This confirms the wobbling phonon excitations in TSD nuclei from 
the microscopic calculations. The result of B{E2) values indicates that the triaxial deformation is 
increasing as a function of spin in the observed TSD bands in ^^''Lu. 

PACS numbers: 21.10.Re, 21.60.Jz, 23.20.Lv, 27.70.+q 



I. INTRODUCTION 



The observation of the wobbling excitations [1, 2] (see 
also Ref. [3, 4]) renewed the interest in the study of 
nuclear rotational motions. The nuclear wobbling mo- 
tion [5] is a quantized motion of the triaxial rotor and 
appears as a multiple-band structure, in which consecu- 
tive rotational bands are connected by strong E2 transi- 
tions with each other. Until now, the multiple rotational 
band structure characteristic to the wobbling phonon ex- 
citation has been observed at the high-spin excited states 
in some Lu isotopes [1-3, 6-10] around ^^'^Lu. Nuclei in 
the Lu and Hf region, including these Lu isotopes, have 
been predicted to be strongly deformed with pronounced 
triaxiality [11-14] at high-spin states, and the associated 
rotational bands are called triaxial superdeformed (TSD) 
bands; i.e., the wobbling structure is composed of these 
TSD bands. In fact the lifetime measurements in these 
nuclei [15, 16] revealed that the rotational E2 transi- 
tion probabilities inside the bands are typically about 
500 Weisskopf units, and those between the bands asso- 
ciated with the wobbling phonon excitations are about 
100 Weisskopf units. These are one of the largest out- 
of-band B{E2) and believed to be the evidence of the 
nuclear wobbling motion. Although several candidates 
of the TSD bands have been observed in even-even Hf 
isotopes [17-21], there is no definite evidence of the wob- 
bling excitation yet. 

The nuclear wobbling motion was first predicted by 
using the simple macroscopic rotor model [5]. The ex- 
perimental data observed in Lu isotopes are investi- 
gated by the particle-rotor model [22-25], because a 113/2 
quasiproton exists in the proton-odd Lu nuclei. It is, 
however, noted that the properties of the observed out- 
of-band E2 transitions suggest that the triaxial defor- 
mation is of the so-called positive-7 shape in the Lund 
convention [26], the sign of which is opposite to that of 
Ref. [5] ; namely the moment of inertia about the shortest 
axis of triaxial deformation is the largest. This conflicts 
with the irrotational moments of inertia which are natu- 



ral for the macroscopic rotor model, where the moment 
of inertia is largest about the intermediate axis. In prin- 
ciple, the triaxiality for the quadrupole shape and that 
for the moments of inertia can be taken differently in the 
macroscopic model; the so-called "7-reversed" moments 
of inertia, i.e., those about the shortest and intermediate 
axes are interchanged, are used in Refs. [22, 23], while the 
rigid-body moments of inertia that are consistent with 
the positive 7 shape are recommended in Refs. [24, 25]. 
Since the moments of inertia are the basic parameters (in- 
puts) of the macroscopic model, microscopic approaches 
are necessary to investigate this problem. 

From the theoretical point of view, how the rotor 
model emerges out of the collective rotational motions of 
constituent nucleons is an interesting and long-standing 
problem (see e.g. Ref. [27]). Most of nuclei are axially- 
deformed in their ground states and the collective rota- 
tion is of one-dimensional nature (rotation about only 
one axis). In contrast, the wobbling excitations corre- 
spond to tilting the main rotation axis from the prin- 
cipal axis of deformation, and the rotation is of three- 
dimensional nature. Therefore, the microscopic study of 
the wobbling motion, which is characteristic in the rotor 
model, is expected to give a new insight to the prob- 
lem of how the individual nucleons form and affect the 
triaxial nuclear rotor. The key to relate the motions of 
nucleons to the collective rotation is the introduction of 
the "body-fixed" or the principal axis (PA) frame [28], 
in which no collective rotations exist. The condition 
of the body-fixed frame defines the constraints to the 
many-body system, and combined with the redundant 
rotor degrees of freedom the original many-body prob- 
lem is recovered. The idea of this "collective coordinate 
method" in quantum many-body systems was put for- 
ward in Ref. [29], especially for the nuclear triaxial ro- 
tor problem in [30, 31]. The rigorous quantum mechan- 
ical treatment of the constraints is a difficult problem, 
but Marshalek showed [32] by taking the small ampli- 
tude fimit of the time-dependent HFB (TDHFB) theory 
that the transformation to the PA frame is straightfor- 
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ward within the order of the random-phase approxima- 
tion (RPA), see also Refs. [33, 34]. The semiclassical 
treatment of the cohective rotational motion based on 
the TDHFB theory was developed in Refs. [35, 36], in 
which the PA frame constraints are taken into account 
properly. The attempts to go beyond the RPA order have 
been done [37, 38] based on the selfconsistent-coUective 
coordinate method [39, 40]. 

In Marshalek's theory [32] , it has been shown that the 
RPA eigenenergy can be written in the same form as the 
wobbling phonon energy in terms of the three moments 
of inertia, which can be derived cither in the classical 
asymmetric top [41] or in its quantum analogy, the triax- 
ial rotor model [5] , by using the small amplitude assump- 
tion consistent to the RPA. It has been also shown [42] 
that the E2 transition probabilities associated with the 
wobbling excitation can be written in the same way as in 
the rotor model [5] if the RPA eigenmodc has a certain 
property. It should, however, be stressed that the actual 
calculations are performed in the cranked mean-field plus 
RPA approach [43-45], which is formulated in the uni- 
formly rotating (UR) frame according to the usual crank- 
ing prescription. It is the introduction of the PA frame 
and the transformation from the UR to PA frame that 
makes it possible to interpret a microscopic RPA eigen- 
modc as the wobbling phonon and to derive the three 
moments of inertia of the rotor microscopically. In the 
theory of Ref. [32] as well as those in [35, 38], the PA 
frame is defined such that the non-diagonal part of the 
(mass) quadrupole tensor should vanish identically [28]. 
It is discussed in Refs. [29, 31] that the PA frame condi- 
tion is a gauge fixing condition and physical observables 
should not depend on its choice; for example, the three 
moments of inertia in the PA (body-fixed) frame them- 
selves are gauge-dependent and are not the physical ob- 
servables. However, in order to obtain a clear physical 
picture and a relation to the macroscopic rotor model, 
the introduction of the specific PA frame condition cor- 
responding to the nuclear rotor is important and useful. 

According to the Marshalek's wobbling theory, micro- 
scopic RPA calculations have been performed for nuclei 
in the Hf and Lu region [46-48], where it was found 
that the RPA eigenmodc that can be interpreted as 
the wobbling phonon exists in each nucleus stably in 
a reasonable range of mean-field deformations (see also 
Ref. [49]). Although it is not confirmed in experiment, 
the possible wobbling excitations in normal deformed nu- 
clei with the negative-7 shapes have been also studied 
previously [42, 50] and more recently [51]. The relation 
between the instability of the wobbling RPA mode and 
the appearance of the tilted axis cranking [52] mean-field 
solution is also discussed [49, 53, 54]. In the previous cal- 
culations [46-48], the calculated out-of-band to in-band 
B{E2) ratio was systematically smaller by about factor 
two than the experimental data, but it was shown in 
Ref. [55] that this is not a problem: The triaxiality pa- 
rameter in the Nilsson potential, which was determined 
by the Nilsson Strutinsky calculation and used in the 



RPA calculation, is too small compared with the value 
estimated by the macroscopic rotor model. If the same 
value as that used in the particle-rotor model analysis is 
employed, the RPA calculation gives correct magnitudes 
of the B{E2) ratio. 

In all the previous realistic calculations of the RPA 

wobbling excitations mentioned above, the Nilsson (mod- 
ified oscillator) potential is used as a nuclear mean-field. 
The problem of the Nilsson potential applied to high- 
spin cranking calculations is well-known: The moment 
of inertia is overestimated due to the potential, see 
e.g. [56]. The simplest way to avoid this problem is to 
apply the Strutinsky renormalization to the angular mo- 
mentum expectation value, which cannot be used in the 
RPA calculation. In fact the absolute value of the mo- 
ment of inertia is about 30% overestimated in the pre- 
vious calculations [46-48]. More fundamental remedy is 
to include a new field-coupling in order to restore the lo- 
cal Galilean invariance broken by the velocity dependent 
potential like the term, which was suggested in [5] and 
formulated in Ref. [57] . Based on this formalism the cor- 
rection to the cranking term was included in the high-spin 
cranking calculation in [58]. Recently, this method has 
been successfully employed in Ref. [59], and the author of 
Ref. [59] claimed that the problem of the Nilsson poten- 
tial is solved. However, when this method is applied, we 
found that no RPA solutions corresponding to the wob- 
bling phonon were obtained for the TSD nuclei in the Lu 
and Hf region. The reason is clear: The method is only 
applied to one component along the cranking axis among 
three possible rotational axes. The field-coupling to re- 
store the local Galilean invariance is realized as a residual 
two-body interaction, which is composed of three compo- 
nents corresponding to the infinitesimal rotations about 
each Cartesian axis, in the original formulation [57] . The 
method to include the correction to the cranking model 
is derived by using the mean-field (HF or HFB) approxi- 
mation to this residual interaction, so that only the com- 
ponent of the cranking (a;-) axis is effective. All the three 
components, especially the perpendicular (y, z-) compo- 
nents, are necessary to be included for the description of 
the wobbling excitation, which is of generic three dimen- 
sional rotation. Thus, the additional residual interaction 
should be included in the RPA wobbling calculation. 

The present work is along the same line as those works 

based on the formalism of the cranked mean-field and the 
RPA, but we use the Woods-Saxon potential as a nuclear 
mean-field to avoid the problem mentioned above because 
there is no velocity dependent part of the potential (ex- 
cept for the spin-orbit term) . Apart from the problem of 
the velocity dependence, the Woods-Saxon potential is 
believed to be more realistic than the Nilsson potential. 
It is desirable to confirm the existence of the wobbling- 
like RPA solutions also in this mean-field potential. This 
is the main subject of the present work. The paper is or- 
ganized as follows: The basic formulation of the wobbling 
RPA theory [32, 42] is reviewed in Sec. II. The theory 
needs to be slightly modified for the general deformed 
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mean-field like the Woods-Saxon potential. The neces- 
sary modification is explained in this section. The results 
of numerical calculations are presented in Sec. Ill, where 
the detailed and systematic analyses are performed. Sec- 
tion IV is devoted to the summary. A part of the present 
work was presented in a conference report [60]. 



II. RPA WOBBLING THEORY 

Although it may be desirable to perform the full self- 
consistent cranked HFB plus RPA calculation by using 
a realistic two-body effective interaction, such as one of 
Skyrme forces or Gogny forces, it is still not very easy 
because all the symmetries except the parity and signa- 
ture are broken in the triaxial superdeformed (TSD) rota- 
tional band. One way is to use a simple mean-field and a 
schematic interaction such as the pairing plus quadrupole 
interaction [27]. We take a different approach: Starting 
from a reliable nuclear mean-field for triaxially deformed 
nuclei, we construct the residual interaction that is suit- 
able to describe the nuclear wobbling motion. In this 
way, we can investigate the properties of the wobbling 
mode for any given mean-fields, e.g., its dependence on 
the triaxiality or the pairing gap. In this section, after 
explaining the Woods-Saxon potential and how to con- 
struct the residual interaction, we discuss the RPA wob- 
bling formalism, which needs to be slightly modified from 
the original one by Marshalek [32]. 



A. Triaxially deformed Woods-Saxon potential 

The triaxially deformed Woods-Saxon potential con- 
sidered in this work is a standard one widely used in the 
study of high-spin states [61-63]: 

Vwsir) = V;(r) + A,o (i^r^) (VK(r)) • (a x ^V). 

(1) 

The first term is the central part and the second is the 
spin-orbit term. Mj-ed = with A being the mass 

number and M being the nucleon mass, and cr is the 
Pauli matrices (the nucleon spin s = ^cr). The explicit 
form of the central potential is 



V 



1 -|- exp(disti;(T")/a) ' 
where the strength V is given by 



V = -VoX [l±K 



N-Z 



-|-proton 
—neutron 



(2) 



(3) 



with Z and N being the proton and neutron num- 
bers, respectively, and distE(T') is the distance between 
a given point r and the nuclear surface S, with a mi- 
nus sign if r is inside S. The nuclear surface S in this 
work is parametrized by three deformation parameters. 



(/32,7,/94), and is defined by the usual radius to solid- 
angle relation, r = 



K=0,±2 



K=0,±2,±4 



(4) 

where Cy{P2,7, Pa) is determined by the volume conserv- 
ing condition, and the coefficients a's are given by 



' 020 = /32C0S7, 

022 = 02-2 = - 
O40 = 5/34(5003^7+ 1) 



71 ^2 Sin 7, 



042 



6 ' 
04-2 



I /34 COS 7 sin 7, 



(5) 



044 = 04-4 = '\J^ Pa sin^ 7. 



The Coulomb potential, which is calculated by assuming 
the uniform charge distribution inside the surface S given 
by the central potential. 



3{Z-l)e^ 



(r 



■dS' 



(6) 



is added to the proton potential. 

The potential is completely specified by the parame- 
ters, Vb, K, Ro, a, and A. In this work we use a set of 
those parameters provided by Ramon Wyss [64], among 
which K and Rq are different in the central and spin- 
orbit potentials. These parameters are determined by 
the requirement that the moments of inertia and the 
quadrupole moments can be nicely reproduced system- 
atically for medium and heavy nuclei throughout the nu- 
clear chart. They are given in Table I. 

It should be mentioned [55] that the triaxiality pa- 
rameter 7 in Eq. (5), for which we denote 7(pot:WS) if 
necessary, is considerably different from the one defined 
in terms of the quadrupole moments, 7(den). For ex- 
ample, 7(den) = 20° corresponds to 7(pot:WS) ~ 30° 
for the case of large deformation such as the TSD bands, 
which gives important consequences on the interpretation 
of the B{E2) observed in the wobbling TSD bands (see 
Ref. [55]). This point will be discussed in more details in 
Sec. HIE. 



B. Residual interaction and RPA in the uniformly 
rotating frame 

The method to construct the residual interaction is 
discussed in relation to the wobbling motion in Ref. [65] . 
The idea is based on the decoupling of the rotational 
Nambu-Goldstone (NG) modes, i.e., the angular momen- 
tum operators Ji {i = x, y, z), within the RPA [66]. The 
same idea is formulated in the context of the particle- 
vibration coupling theory [5]. Although the following dis- 
cussions arc essentially the same as those of Refs. [32, 33], 
we recapitulate them because some of their results are 
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TABLE I: The parameters of the Woods-Saxon potential used in this work [64]. Other physical constants used are 
e^/Rc=137.03602, he = 197.32891 MeV-fm, and Mc^ = 938.9059 MeV. A in the table denotes the mass number of nucleus. 



Vo (MeV) 






Roc (fm) 


-Roso 


a (fm) 


Aso 


53.7 


0.63 


0.25461 


1.193^1/^ + 0.25 


0.969 X iloc 


0.68 


26.847 



used in the later discussions in order to make clear the 
point that should be modified from the original RPA wob- 
bling theory [32]. It should be noted that the signature, 
the TT rotation about the cranking axis (x-axis), is a good 
quantum number, and only the part of the RPA equa- 
tion which transfer the signature quantum number a by 
Aa = 1, the so-called signature (— )-part, is relevant to 
the wobbling motion. Therefore, we consider only the 
corresponding part of the residual interaction. 

Starting from the general mean-field hamiltonian, h, 
which, in the present work, is composed of the kinetic 
energy, the Woods-Saxon potential, the Coulomb poten- 
tial, and the pairing part (see Sec. Ill A), the relevant 
part of the residual interaction is given by 



Hres = -lxyF^-lx^{^F,f, 

where the operators Fy and are defined by 
Fy = [h, iJy] , iF^ = i[h,J;,]. 



(7) 



(8) 



The strengths of the interaction are determined by the 
decoupling condition of the NG modes, Jy and Jz, for 
the total hamiltonian, H = h + -ffrcs, 

[H, Jy,z]RPA = [h + i?res; Jy,z]RPA = 0, (9) 

where the subscript RPA means that the commutator is 
evaluated within the RPA order, leading to 



VXy = {{[h,iJy],iJy]}, 

l/Xz = -([[/l,Jz],Jz]). 



(10) 



Here (and hereafter if not stated explicitly) the expecta- 
tion values are taken with respect to the cranked mean- 
field state (the generahzed product state) |<I>(ci;rot)), 
which is a vacuum state of the quasiparticles in the uni- 
formly rotating frame about the x-axis. 



(11) 



where Wrot is the rotational (cranking) frequency, Ea is 
the energy of the quasiparticle eigenstate \a) in the rotat- 
ing frame (quasiparticle routhian), and a|j is its quasipar- 
ticle creation operator. It should be stressed that the vac- 
uum mean-field state, |$(a;i.oi))! determines the strengths 
in contrast to the conventional case, where the two-body 
interaction, such as the QQ force, is given and it deter- 
mines the vacuum mean-field state. This is very impor- 
tant to guarantee the NG mode decoupling: It is usually 



a numerically demanding task to satisiy the selfconsis- 
tency between the interaction and the vacuum mean-field 
state in such a level as to exactly satisfy the NG mode 
decoupling. 

The creation operator of the RPA eigenmode excited 

on the cranked vacuum state |<I>(ciJrot)), i-c., the one in 
the uniformly rotating (UR) frame, is given as a linear 
combination of the two-quasiparticle excited states, 

= [^niaP)aial + <pn{a/3)ai3aa , (12) 

a</3 

and is obtained by solving the RPA equation in the UR 
frame, 

[H\jK, ^^]rpa = w„X^, Hub. = H — uJrotJx, (13) 

where w„ is the corresponding RPA eigenfrequency. Here 
and hereafter we use h = 1 unit for simplicity. 

The operators which play a key role for defining the 
principal axis (PA) frame discussed in the next subsec- 
tion, are the non-diagonal part of the (mass) quadrupole 
tensor operators. What is relevant to the RPA wobbling 
theory is their signature (— ) part, 



Qy — Q21 ^ ~ 



Qz — Q22 ^ — * 




(14) 



0=1 



where the 2;-axis is chosen as a quantization axis of the 
signature classified quadrupole tensor. 



^2K 



1 



V2(l + Sko) 



{Q2K±Q2-k), (K = 0,1,2). 

(15) 

They are important also because the out-of-band E2 
transition probabilities are calculated by their RPA ma- 
trix elements. The part composed of the two quasipar- 
ticle excitations only contributes within the RPA order, 
and it is denoted with the superscript {A) in the foUow- 
ings; 



' (Oy)^^^ = XI «!/("/5)(al«^ + a/3aa), 



a</3 



(<9J^^^ = X 9^("/3)(44 ~ "'0^°')^ 

a<l3 



(16) 



where wc used the phase convention [67] that the two- 
quasiparticle matrix elements in this and the following 
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Eq. (17) are all real, qy^z{a(3) = q* .,{al3), jy,z{ct(3) = 
j*_^(a/3). As it can be easily checked, the signature (— ) 
part of the angular momentum operators, 



a</3 
a</3 



(17) 



combines into the RPA eigenmode corresponding to the 
NG mode in Eq. (13), 



^n=NG 



1 



{zJy)^^^ + {JzY^^], (18) 



with the eigenfrequency a;„=NG = ujrot- Here the expec- 
tation value of the angular momentum along the cranking 
axis can be written as 



a<l3 



(19) 



In the same way, for the operators Fy z in Eq. (8), 

a<l3 

(J^,)(^) = ^/,(a/3)(at4-a^a„), 

a</3 

and by the definition their matrix elements are expressed 

as 



(20) 



fy{al3) = Eafjjy{al3) ~ uj^odziafS), 



(21) 



with Eai3 = Ea + E/j, and the force strengths as 

ct<0 a<0 

Xz' = 2 J2 Ea0fz{ap) - 2a;,ot ^ jy{aP)jz{ap). 

(22) 

Since the residual interaction is of the multi- 
component separable type, it is straightforward to solve 
the RPA equation by means of the dispersion matrix 
technique [44, 45]. The forward and backward ampli- 
tudes in Eq. (12) are given by 



Ea0 — 

-fy{aP)Xy^yin) + fz{aP)Xz^z{n) 



(23) 

where the RPA transition matrix elements for the oper- 
ator Fy_z arc defined by 

' J-^(n) = {[Xr,,Fy\) = ^(V'n(a/3) - K{ap))fy{ap), 

a<l3 

Tz{n) = ([X„,F,]) = ^ (V'„ (a/3) +</.„ (a/3)) /,(«/?), 

a</3 

(24) 



which are obtained by solving the following linear equa- 
tion 



RyyioJn) — Xy^i Ryzi^n) 

Rzy{uj„), RzzioJn) - Xz^ 



XyJ'vin) 
Xz^zin) 



Here the quantities in the matrix are given by 

RyyiLo)^D(^^fy,fy;cu), R z z (co) = D^^^ {f z , f z', uj) 
Ryziu) = Rzy{u) = D(^\fyJz;Lo), 

where the response functions D^^'"^^ are defined by 

2E^^f{a/3)g{ap) 



= 0. 

(25) 



(26) 



a<l3 



D('\f,g;uj)^J2 



a<0 



p2 , ,2 

2ujf{af3)g{a(3) 



(27) 



Thus the RPA eigenfrequency a; = a;„ is obtained by 
solving the dispersion equation 



det 



-R^/!/M-X^/^ Ryz{i^) 
Rzy{oj), Rzz{uj)-Xz^ 



(28) 



and the norm of the corresponding solution {Xi^iin), i = 
y, z) is determined by the normalization condition 

Eo</3(V'SM)-0^M)) = 1- 

In the case of the residual interaction (7) with (8), the 
NG mode (18) can be explicitly separated out from the 
equation. By inserting Eqs. (21) and (22) into Eq. (25), 
the matrix in Eq. (25) can be written as 



UlnBy{cOn) — Wrot -^zi^^n), W„ (a;„ ) — Ul^ot By{uin) 
U>nAy{u)n) - LJrotBzii^n), ^nBz{^n) " <^Tot-^y{^n) 



with the definitions, 

Ay(uj) = W,.ot(jx - Jy(w)) + Ujjyziuj), 
Az{uj) = UJrot{Jx - Jzi^j)) + U)Jyz{oj), 
By{uj) = Ujjy{uj) - U>TotJyz{'^), 
. Bz{l^) = t^Jz{i^) - t^TotJyzii^), 



(29) 



(30) 



and 



X = Jy{co) = D^'\jy,3y;cj), 

Jzico) = D(i)(j„j,;a;), Jyzico) = D(^\jy,jz;io). 



Then Eq. (25) can be cast into the form 



0, Wrot + 

Az{uJn) + By{uJn), Bz{uJn) + Ay{u)„) 
Az{uJn) - By{0Jn), Bz{0Jn) - Ay{uJn) 



(31) 



Xy^yin) 
Xz^z{n) 



0, 



(32) 



6 



from which the a;„=NG = '-^rot NG solution is apparent, 
because by using Eq. (18) the NG amphtudes {xi^i, i = 
y, z) reduce to 



XyJ'yin = NG) 
XzJ^z(n = NG) 



1 



Thus, for the non-NG solutions, ojn 7^ ijJiot 



(33) 



Xy^y{n) 



0, 



(34) 



and the dispersion equation for them is 



det 



By{u),Ay{u) 



0, 



(35) 



which was first derived in Ref. [33]. Alternatively, 
Eq. (25) can be cast into an another form 



Ay{uJn),Bz{uJn) 
By{u}n),Az{u)n) 



0, 



with the definitions, 



fiy{n) ^ LUnXy^y{n) - i0^otXz^z{n) , 
flzin) = (jJnXz^zin) - UJrotXyJ^yin)- 



(36) 



(37) 



Again, this equation contains the trivial NG solu- 
tion (33), i.e., fiy(n) = f2z(n) — 0, with a;„=NG = 
Wrot- For the non-NG solutions, fiy{n) / flz(n) = 
-Bz{ujn)/Ay{LOn) = -A^(a;„)/By (w„), therefore, by 
defining 



Bz{cOn)' 
Az{u)n) 

Eq. (36) can be further changed into 



jf{n) = Jy{Wn)-Jyz{u^n) 



J fin) = Jz{uJn) - Jyz{^n) 



(38) 



^.ot{Jx- Jf{n)), LOnJfin) 
WnJfin), w,ot(Jx - Jf{n)) 



{2y{n) 
Qz{n) 



0, (39) 



from which the wobbling phonon energy is formally 
solved, 



Jf{n)Jf{n) 



(40) 



This equation is equivalent to Eq. (35) and was first de- 
rived by Marshalck [32]. 

In the original formulation of Ref. [32] , the quantities 
Qy^zin) and jfzi''^) ^^"6 interpreted as the amplitudes 

of y, z-componcnts of the angular frequency vector and 
of moments of inertia, respectively, in the principal axis 



(PA) frame, when the n-th RPA mode is excited. How- 
ever, this is not the case for the general deformed mean- 
field like in the present work, which is the subject of the 
next subsection. 

For the sake of completeness, the E2 and Ml transi- 
tion probabilities from the n-th RPA phonon excited ro- 
tational band to the vacuum band are calculated, within 
the lowest order in the 1// expansion consistent to the 
RPA order [43], by 

B{E2- In^I± Ivac) « | ([Of±i, Xt])|^ (41) 
B(Ml;/n^7±lvac)« \{[^H±l, Xl])\^ . (42) 

Here the electric E2 and Ml operators defined with re- 
spect to the a;-axis, and {u = ±1) are given by 



I 



i2l 



Ml±l 



(43) 



(44) 



where the electric quadrupole operators, Q2K (^^ = 
1, 2) and Q^^, arc defined in the same way as in Eqs. (15) 
and (14), but only the proton contributions are included 
and the electric charge e is multiplied. The magnetic 
moment operator /x^ (A; = x, y, z) is given as usual. 



IJ'k = y ^IJ'N^{gl^^lk+ 9i''^Sk)a (r = 7r,i/). (45) 

a=l 

As for the in-band stretched E2 transition probabilities, 

those in the phonon excited band and the vacuum band 
are the same in the RPA order and are calculated by 



B(£2;7^7±2)« i|(Qf±2)r 



(46) 



C. Interpretation of the RPA eigenmode in the 
principal axis frame 

All the calculations of the observable quantities such as 
excitation energies and electromagnetic transition prob- 
abilities can be done within the RPA in the UR frame as 
has been reviewed in the last subsection. It is, however, 
necessary to go into the PA frame in order to interpret 
the obtained RPA eigenmode and compare its property 
with that of the wobbling motion, which is introduced in 
the macroscopic rotor model. Furthermore, the transfor- 
mation is inevitable to investigate how much the angular 
momentum vector, or the angular frequency vector, tilts 
or wobbles around the main rotation axis. This has been 
done by Marshalck [32] , and reviewed from a slightly gen- 
eral view point in Ref. [42] in the light of the quantum 
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theory of collective coordinates [29, 31]. In this subsec- 
tion, we recapitulate some results from Ref. [42] , in order 
to show that a slight modification of the original for- 
mulation in [32] is necessary for the general deformed 
mean-field. 

The transformation from the laboratory frame to the 
UR frame is an unitary transformation so that the RPA 
in the UR frame satisfies all the usual properties [27], 
e.g., the completeness relation. Therefore, the two- 
quasiparticle excitation part of the quadrupole operators 
in Eq. (16) can be expanded in terms of the RPA eigen- 
modes, 



mail 



(47) 



n:all 



where the amplitudes Qy,z{n) are calculated by 

' Qy{n) = {[X^,Qy]) = ^(t/'„(a/3) -</.„(a/3))gj,(a/3), 

a</3 
a<l3 

(48) 

In particular, the NG mode contribution is contained in 
Eq. (47), 



Q^(n = NG) = 



Q,{n = NG) = 



2a^i?2 



(49) 



^/2(X) 



where i?^ = {Y.^^^{r'^)a) is the mean square radius and 
ay{az) describes the deformation around the y (2;)-axis, 



2ayR'^J^£{Y,{x^-zX) 



o=l 



= {[iJy, Qv\) = 2 3v{.ali)qy{ap), 



a<l3 



(50) 



a—1 



In contrast, the two-quasiparticle excitation part of the 
angular momentum operators are composed purely of the 
NG mode, 



(51) 



X- 



NG 



ATng)- 



the RPA wobbling mode is excited, from which the mo- 
ments of inertia are naturally introduced. Wo do not go 
into the full details of the theory, but briefly mention 
the procedure. As it is well-known, the (three dimen- 
sional) cranking prescription can be viewed as a result of 
the constraints, Ji = li {i = x,y,z), where li are time- 
dependent constants in the semiclassical theory [35, 38], 
or are the total angular momentum operators in the rotor 
model (i.e., the differential operators with respect to the 
Euler angles) in the quantum theory of collective coor- 
dinates [30]. This kind of constraints are called the first 
class constraints in the general theory [68], and one has 
to impose the second class constraints, i.e., the PA frame 
condition in this case, in order to develop the consistent 
framework of the quantum [29] or the classical canoni- 
cal [40] theory. The second class constraints can be, in 
principle, arbitrarily chosen (the gauge condition), and 
the physical observables should not depend on the choice. 
However, the actual choice is important for the physical 
picture; recall, e.g., the choice of the center of mass co- 
ordinate in the case of the translational motion. Here we 
take the following PA frame condition [28, 32, 35]; the 
non-diagonal part of the quadrupole tensor should vanish 
identically: 



(Qy)pA = 0, (Q,)PA = 0, 



(52) 



(the constraint for the .x-component, i.e., {Qx)pa = 

{Q2i^)pA = 0, is automatically satisfied in the RPA or- 
der), where the subscript PA denotes that the quantity 

is expressed in the PA frame. 

This condition relates the collective coordinates, i.e., 
the Euler angles in the present case, to the nucleon de- 
grees of freedom. Although it is a difhcult task in general, 
it has been shown in Ref. [32] that this relation can be 
solved within the small amplitude limit (RPA). The UR 
frame is transformed to the PA frame by the rotation in 
terms of the dynamical Euler angle variables ©, 



(o)PA = c/(e)0[/~i(e), 



(53) 



for an arbitrary operator O. Here we omit the subscript 
UR for the quantity in the UR frame. Within the RPA 
order [32], 



U{Q) w 1 — i6[smtjjjy + costjjJz), 



(54) 



where O = (tpjOjCp) are Euler angles with respect to the 
X-axis (sec Ref. [32]), and the small amplitude limit, <C 
1, is used. From Eqs. (54) and (53), 

(0)PA ^ O + 0smij{[O,iJy]) + i9cos^{[0,J,]), (55) 

and applying this to Eq. (52), one obtains 



2a^R? 



, iO cos i}) — 



(56) 



Now let us consider the behavior of the angular mo- 
mentum vector and the angular frequency vector when 



In this way the collective Euler-angle variables can be ex- 
pressed directly by the nucleon degrees of freedom within 
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the RPA order, i.e., the constraint conditions are approx- 
imately solved (the remaining angle </> does not play any 
role within the RPA order). Applying the transforma- 
tion (55) to the angular momentum operators, 



{A) 



(57) 



^PA 



J: 



{J,} r,(n)(xt+X„), 



namely, their amplitudes associated with the n-th RPA 
eigenmode in the PA frame are 



f {Jy)pAin) = {[Xn, {iJy)pA]) = -{Jx)rz{n) 
\ {Jz)pA{n) = {[X„,{Jz)pa\) = {Jx)ry{n), 

with the definitions. 



(58) 



which describe the shape fluctuations (ratios of the dy- 
namic fluctuations amplitude to the static deformations). 
Thus, the contribution of the NG mode, which will be 
identified as the collective coordinate, is driven out and 
only the intrinsic degrees of freedom remain in the PA 
frame. From Eq. (57) it is clear that the angular momen- 
tum vector fluctuates or wobbles about the main rotation 
axis (x-axis) , when the n-th RPA phonon is excited. The 
quantity rj,_z(n) represent their fluctuations; 



. . _ (J.)pA(n) _ (Jy)pA(n) 



(60) 



Therefore these quantities should be small Vy^z {n) <C 1 in 
order for the RPA to be valid. Note that all the ampli- 
tudes Qy^z{n) and ry^z{n) are real because of the present 
phase convention. 

The angular frequency vector operators (i7j, i = x,y, z) 
in the PA frame are introduced as Lagrange multipliers 
of the first class constraints, Ji = U, and the dynamical 
time dependence of nucleon operators in the PA frame is 
generated by the hamiltonian in the PA frame. 



Hp A — H — f^x'^x — ^y^y — ^zJz 
~ -ffuR — ^yJy — ^zJz, 



(61) 



within the RPA order, because f}^ ~ t^rot- The angular 
frequency variables {^y,z) are determined by the consis- 
tency condition of the PA frame (gauge) condition, i.e., 
Eq. (52) should hold in arbitrary time. 



IJ^{QV,.)PA = ([Qy,.,i?PA])pA = 0. (62) 



Here {[Qy, HpA])pp^ ~ {[Qy, HuB.])p^ + if2y{[Qy,iJy]), 
( [Qz , Hpa] ) pa ~ ( [Q. , H^b])pa - < IQ- ' "^J > ' and 
HvR ~ J2n:a.n '^niXlXn + 5) within the RPA order, and 
using the transformation (55) for {[Qy,z, -f^UR])pA) '^^ ob- 
tain 

{iny)^^^=- {tOnry{n)+u;rotrz{n)){Xl-Xn), 

(I2,)(^)= Y {i^nr^{n)+u,otry{n)){xl + X^), 

(63) 

i.e., their amplitudes associated with the n-th RPA eigen- 
mode are 

f ■^yW = {[Xn,i^y]) = -{i^nryin) + u;rotrzin)) , 
1 ^z{n) = {[X„, Oz]) = i^nrz{n) + (JrotTyin). 

(64) 

Again, the NG mode contribution cancels out completely. 
From Eq. (63) the angular frequency vector also fluctu- 
ates about the main rotation axis, and the moments of 
inertia in the PA frame are naturally introduced; 



Jf{n) 



(Jy)pA(n) 

Oyin) 

{Jz)pA{n) 



(65) 



f2z (n) Wnrz (n) + uj^otTy (n) 



It is clear that the moments of inertia thus defined are 
intimately connected to the PA frame condition (52) and 
determined by the RPA transition amplitudes of the Qy^z 
operators in Eq. (48). 

Prom Eq. (65) the RPA eigenenergy and the ^^^^(n) 
amplitudes can be represented in terms of these moments 
of inertia [32, 42]; 



■ Wrot 1 



\Jx-Jf{n)){Jx-Jf{n)) 
Jf{n)Jf{n) 



(66) 



i.e., the well-known formula in the rotor model [5] is re- 
covered, and 



ry{n) 
rz{n) 



1 

'"y2(X) 
1 



l/Jf{n) - 1/Jx 

l/jfin) - l/J^ 
llJf{n) - l/J^ 



-,1/4 



1/4 



(67) 



where the quantity c„ and the sign cr„ are defined by 



c„ = [sign of ry{n)\ x ^2{Jx)\ry{n)rz{n)\, 
an = sign of {ry{n)rz{n)). 



(68) 



Combining Eqs. (67), (59) and (41) with Qf^ being re- 
placed to {eZ/A)Q'2v, the rotor model expression for the 

out-of-band B{E2) is also recovered, if (7„ = + and 
cl = 1 are satisfied for the excitation of the n-th RPA 
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cigcnmode. This was first shown in Ref. [42], and the con- 
dition, iT„c^ = 2{Jx) ry{n)rz{n) = +1, is used to identify 
the wobbling mode in a recent pubhcation [51]: Note that 
(7„ = — and c„ = — 1 for the NG mode. 

It must be emphasized that the discussion above is 
general and does not depend on specific choices either 
of residual interactions or of nuclear mean- fields. In 
the original formulation of Marshalek [32] , the arbitrary 
spherical mean-field and the QQ (plus the monopole pair- 
ing) force are assumed. It can be easily checked that the 
formulation is still valid for the selfconsistent anisotropic 
harmonic oscillator potential model [5] , for which the in- 
duced residual interaction is the doubly-stretched Q"Q" 
force [69-72]. In such cases, the deformation is purely 
of quadrupole type, and the selfconsistent mean-field h 
leads for Fy^^ in Eq. (8) to 

and the quantities i7y_^(n) and Jy%Xn) introduced in 
Eqs. (37) and (38) coincide with f2y^z{n) and Jy^^i^) i'^ 
Eqs. (64) and (65), respectively. However, this is not the 
case generally, and the moments of inertia should not 
be calculated by Eq. (38) but by Eq. (65) for the gen- 
eral mean-field: This is the point, in which the original 
formulation of Marshalek should be modified. As is dis- 
cussed above, the moments of inertia are determined by 
the PA frame condition. Therefore, a possible alterna- 
tive choice of the PA frame condition is {Fy^z)pA = 
in place of (Qj,.z)pa = 0. Then the original theory is 
recovered, but now the rotor model expression for the 
out-of-band B{E2) cannot be derived, because the oper- 
ators Fy^z are not related strictly to the quadrupole op- 
erators Qy^z for the general mean-field. We use Eq. (52) 
as the PA frame condition for the realistic calculations 
presented in Sec. III. 



D. Remarks on the NG mode decoupling 

In the realistic calculations, it is not an easy task to ex- 
actly realize the NG mode decoupling, even though the 
microscopic hamiltonian satisfies the decoupling condi- 
tion (9). For example, one has to use a truncation of the 
model space for diagonalization or a finite mesh size to 
solve differential equations of the single-particle states, 
which inevitably breaks the condition. When starting 
from some effective two-body force and the mean-field is 
determined by the HF (or HFB) procedure, a very accu- 
rate achievement of the selfconsistency between the in- 
teraction and the resultant mean-field is required. Prac- 
tically, it often happens that even the decoupling condi- 
tion (9) does not rigorously meet. For example, if one 
takes the Nilsson potential and the QQ force (or the 
doubly-stretched Q"Q" force) as the total hamiltonian, 
then the rotational NG modes do not decouple exactly, 
because the general Nilsson potential [26] contains the 



"hexadecapole" deformed potential (the 64 deformation) , 
whose effects cannot be dealt with the simple QQ (or 
Q"Q") force. Even if the £4 deformation is neglected, 
the / • s and terms are those defined in the stretched 
coordinate in the standard Nilsson potential, so that they 
are not strictly spherically invariant. 

The previous calculations of the TSD wobbling exci- 
tations in Refs. [46-48] suffer from this problem; even 
though the £4 deformation was neglected in them, the 
stretched I ■ s and terms are used. Then Jy^zi^) 
Eq. (38) does not coincide with Jy^zi^) i^i Eq. (65) rig- 
orously. However, we have checked that the effect of the 
model space truncation is much severe, and the effect of 
the stretched I ■ s and terms is very small; jfy5('^) 
and J^^zin) 

agree within 1% if enough numbers of the 
oscillator shells are included. In the recent calculations 
in Refs. [51, 59], the non-stretched I ■ s and terms are 
used in the Nilsson potential without the £4 deformation 
consistently, but the additional correction term, which 
recovers the Galilean invariance broken by the velocity 
dependent part, breaks the spherical symmetry. There- 
fore, in all these calculations, even if £4 — 0, the NG mode 
decoupling was not strictly realized. In Refs. [51, 59], one 
of the strengths of the Q"Q" force is adjusted so that the 
RPA spectra contain the eigenenergy ujn — t^rot, but then 
the corresponding eigenmode solution might not gener- 
ally be the correct NG mode, Eq. (18), although it is 
claimed that such effect is small [51, 59]. 

If the monopole (or seniority) pairing correlation is in- 
cluded in the mean-field, it can be an another source 
that disturbs the actual NG mode decoupling. Since 
the monopole (or seniority) pairing interaction gives di- 
vergent results in the full model space, some truncation 
scheme is necessary. Usually the energy truncation of the 
single-particle states is used, in which not all of the mag- 
netic substates (m) for the orbital with a given angular 
momentum (7) arc included in the pairing model space, 
and then the rotational invariance is broken. We have 
checked that this effect is very small if the pairing model 
space is not too small. 

The method to construct the residual interaction that 
restores the rotational symmetry in Sec. II B can be used 
to remedy (or avoid) this problem generally. Since the 
wobbling motion is intimately related to the rotational 
NG mode, we believe that it is important to respect its 
decoupling, and the method in Sec. II B is a simple and 
flexible way to realize it. 



III. RESULTS OF NUMERICAL 
CALCULATIONS 

The purpose of the present work is to study the wob- 
bling excitations of the triaxial superdeformed (TSD) 
bands in the Lu and Hf region. We first discuss the re- 
sult of the mean-field calculations based on the new pa- 
rameterization of the Woods-Saxon potential explained 
in Sec. II A, and then the result of the microscopic RPA 
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calculations. 



A. Mean-field calculations 

Our mean-field hamiltonian is 

hr = tr + VwsAr) + Sr,nVcoul{r)-Ar{Pl+Pr)-XrNr, 

where t = tt, distinguish the proton or neutron part, t 
is the kinetic energy term, Vws,t('") is the Woods-Saxon 
potential in Eq. (1), Vcoui('") is the Coulomb potential 
in Eq. (6), and A are the monopole pair transfer op- 
erator and pairing gap, and N and A are the number 
operator and the chemical potential, respectively. We 
have developed our own code, based on the program 
of Ref. [73], for the cranked Woods-Saxon calculation 
with pairing correlations included. The diagonalization 
is performed in two steps; first the Woods-Saxon (and 
Coulomb) potential {A = and Wrot = 0) is diagonalized 
in the anisotropic harmonic oscillator basis, and then the 
cranked pairing problem (the HFB equation) is solved in 
the Woods-Saxon basis obtained in the first step. 

The matrix elements of the Woods-Saxon (and 
Coulomb) potential with respect to the anisotropic oscil- 
lator basis are specified by the oscillator quantum num- 
bers, {nx,ny,nz). They are evaluated by using the re- 
currence relations of the Hermite polynomials [73] start- 
ing from the diagonal and subdiagonal matrix elements, 
which are calculated by the Gauss-Hermite quadrature 
formula. The three oscillator frequencies, {ujx,'jjy,ujz)i 
arc chosen to be inversely proportional to the mean radii 
of the X, y, z-directions, which arc calculated geometri- 
cally by using the uniform density distribution based 
on the shape specified in Eq. (4) for given deforma- 
tion parameters (/32,7,/34). The magnitude of the fre- 
quencies, {iJxiOyiJz)^/^ , is chosen to be 1.25 x loq, with 
^0 = 41/^1/3 MeV. We use the truncation of the oscil- 
lator shells, iVmax = {nx+ny+nz)ma.-yi = 12, in most of the 
calculations presented in this work. As for the Coulomb 
potential (6), it is evaluated by the Gauss- Legendre nu- 
merical integration if viimQR{Q) < r < ma,XQR{n), and 
by using the multipole expansion of \r — r'\^^ up to 
^max = 16 if otherwise. Since we use the monopole pair- 
ing interaction, the truncation of the pairing model space 
is necessary. Woods-Saxon single-particle orbitals whose 
energy satisfies [ej — A| < fuvo are included in the pair- 
ing model space. 

The main reason why we use the Woods-Saxon poten- 
tial rather than the Nilsson potential is that the crank- 
ing moment of inertia J^x in Eq. (31) is overestimated 
unless the Strutinsky renormalization [26] is not per- 
formed for the expectation value (J^). In order to show 
how this overcstimation is greatly improved, we com- 
pare the cranking moment of inertia calculated by us- 
ing the Woods-Saxon and the Nilsson potentials in Ta- 
ble II. The cranking moment of inertia is sensitive to the 
shell structure and the pairing correlation, so that we 



tabulated the Strutinsky smoothed moment of inertia, 

i.e., J'x = (>/a;)/wrot) without pairing for selected nuclei 
and deformations in the rare earth region. This iner- 
tia J^x is Wrot-dependent but its dependence is known to 
be weak [56] (we have checked that the variations are 
less than 2% in the present examples), so that those at 
t^rot = 0.6 MeV/fi. are tabulated. We also include the cor- 
responding rigid-body values calculated by assuming the 
uniform density distribution with the radius 1.2 x A^^^ 
fm. For the calculations of this J'x and of the selfcon- 
sistent deformation of the TSD minima discussed in the 
followings, we employ the cranked Nilsson Strutinsky cal- 
culation of Ref. [74], where the ANosc = ±2 coupling 
arising from the cranking term is included in an approx- 
imate way. The parameters of the present Woods-Saxon 
potential are given in Table I, and those of the Nilsson 
potential are taken from Ref. [75]. The parametrizations 
of shape in the two potentials are different. We convert 
those of the Nilsson potential (£2,7,64) [75] to those of 
the Woods-Saxon potential (/32,7, /34) by requiring the 
same mean square radii of the x, y, z-directions, which 
are calculated geometrically by the uniform density dis- 
tribution. We use this conversion between (,52, 7, A) and 
(^2, 7, £4) throughout this work. As it is clear in Table II, 
the overestimation with the Nilsson potential by 25—27% 
is much more reduced to 6—7%; the remaining deviations 
are mainly due to the spin-orbit potential. 

TABLE II: The Strutinsky smoothed cranking moment of in- 
ertia in unit of ?),^/McV for various nuclei and deformations. 
The corresponding rigid-body values are also included. The 
deformation parameters (/32,7,/34) are those used in Eq. (5). 
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The result of Table II indicates that the single-particle 
orbitals of the present Woods-Saxon potential have cor- 
rect radial size, which is also apparent from the values of 
the radius parameters in Table I. In fact we have checked 
the proton and neutron density distributions calculated 
by using the present Woods-Saxon potential are quite 
similar to those obtained by the Hartree-Fock calcula- 
tions using the standard Skyrme and Gogny forces. 

In order to study the wobbling excitations of the TSD 
bands in the Lu and Hf region, we have to choose suit- 
able deformation of the mean-field. The standard way 
is to search the minima in the Strutinsky method, but 
we are not yet able to perform reliable cranked Woods- 
Saxon Strutinsky calculations of the potential energy sur- 
face with fixed spin values. The main purpose of the 
present work is to confirm the existence of the wob- 
bling excitation modes in the microscopic Woods-Saxon 
RPA calculation. Therefore, we rely on the results of 
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FIG. 1: The single-particle energies as functions of /32(/34) and 7. The upper panel is for the proton energies and the lower 
panel for the neutron energies. In each figure the left part is drawn by changing (52 {Pi) from to 0.42 (from to 0.02) with 
keeping 7 = 0, while the right part by changing 7 from 0° to 60° with keeping {P2, pi) = (0.42, 0.02). The solid (dashed) lines 
are used for positive (negative) parity orbitals. 



the cranked Nilsson Strutinsky calculations for choosing 
the deformation parameters. The obtained deformation 
parameters for the yrast TSD band of ^^"^Lu, i.e.. the 
band with the parity and signature (7r,a) = (+,1/2), 
change very little as functions of spin; they are about 
£2 = 0.41 - 0.38, 7 = 19° - 21°, and 64 = 0.03 - 0.05 



in the spin range r = 25/2+ - 97/2+. Taking the 
(£2,7,64) = (0.39, 20°, 0.04) as typical ones, the con- 
verted values are {P2,1,M = (0.42, 18°, 0.02). We 
choose these values as reference values in the following 
calculations. 

Now we discuss the single-particle properties of the 
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present Woods-Saxon potential. In Fig. 1 the depen- 
dences of the proton and neutron single-particle energies 
on the deformation parameters (32(^4) and 7 are shown. 
It is clearly seen in the present Woods-Saxon potential 
that sizable shell gaps at Z = 70 ~ 74 for protons and 
A'' = 92 ~ 96 for neutrons are existing at the triaxial de- 
formation 7 w 20°. When cranked, these shell gaps are 
responsible for the TSD bands systematically observed 
in nuclei in the Lu and Hf region [13]. Because of the 
relatively large deformation, P2 > 0.4, the iVosc = 6 pro- 
ton orbitals resulting from the 7rii3/2 state come down 
near to the Fermi energy of the Z w 71 system, the low- 
est state of which is occupied in the yrast TSD band in 
Lu nuclei. It is well-known that the occupation of this 
down-sloping and highly- alignable orbitals increase the 
deformation from the normal deformed, (32 ~ 0.2 and 
7^0°, to the triaxially superdeformed shape, /32 ~ 0.4 
and 7 « 20°. 

PROTON 




0.2 0.3 0.4 0.5 0.6 
NEUTRON 




0.2 0.3 0.4 0.5 0.6 



FIG. 2: The quasiparticle routhians as functions of the ro- 
tational frequency suitable for ^^''Lu. The fixed deformation 
parameters {132, 'y, Pi) = (0.42, 18° , 0.02) and the fixed pair- 
ing gaps = Au = 0.5 MeV are used. The upper panel 
is for the quasiproton energies and the lower panel for the 
quasineutron energies. The solid (deished) lines are used for 
signature a = 1/2 (a = —1/2) orbitals. 

The proton and neutron pairing gaps, and A^,, re- 
spectively, are also important mean-field parameters, for 
which we know very little in the TSD bands. There- 



fore, we use constant values for qualitative investigation 
of the TSD bands and of the wobbling excitations on 
them. We will discuss the result where the pairing gaps 
are changed in the following Sec. HID. The chemical 
potentials and are always adjusted to reproduce 
the correct proton and neutron numbers, {Nt^) = Z 
and (iVy) = N. We show the quasiparticle energies 
(routhians) obtained by diagonalizing the cranked mean- 
field hamiltonian, Eq. (11), in Fig. 2, as functions of 
the rotational frequency. Constant pairing gaps, = 
A^ = 0.5 MeV, are used with the reference deformation, 
(/32, 7, A) = (0.42, 18°, 0.02) suitable for ^^^Lu. It is clear 
that there is no crossing of the neutron quasiparticle be- 
fore Wrot ~ 0.6 MeV/h. In contrast, one highly- alignable 
a = 1/2 proton quasiparticle is clearly seen; its main 
component is the A^osc = 6 oscillator shell coming down 
from the 77^13/2 orbitals. This A^osc = 6 quasiproton is al- 
ways occupied in the lowest TSD configuration of ^^^Lu 
and contributes to the alignment of the angular momen- 
tum about the cranking axis. After lu\-oi ~ 0.6 MeV/fi, 
there are many quasiparticles, both proton and neutron 
ones, crossing the zero energy, which changes the configu- 
rations of the yrast and near yrast states and also affects 
the result of the RPA calculations. 

In the following we will show the result of calculations 
mainly for ^^'^Lu, for which most complete experimental 
data are available. 
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FIG. 3: The cranking moment of inertia in Eq. (31) for 
^®^Lu as a function of the rotational frequency. The fixed 
deformation parameters (/32, 7, (^4) = (0.42, 18°, 0.02) and the 
fixed pairing gaps At, = Av = 0.5 MeV are used. Experimen- 
tal data are taken from Ref. [6]. 



Next we compare the calculated and measured crank- 
ing moment of inertia for the yrast TSD band of -'^^^Lu 
in Fig. 3, and those of in-band B{E2) in Fig. 4, respec- 
tively, as functions of the rotational frequency. Here we 
use the standard procedure [76] {K = 1/2 is assumed to 
transform the total angular momentum quantum num- 
ber to the aligned angular momentum) for obtaining the 
measured cranking moment of inertia. For the calcula- 
tion in these figures all the mean-field parameters are 
assumed to be constants, (/32,7,/?4) = (0.42, 18°, 0.02) 
and Z\7r = Aij = 0.5 MeV, as a typical example. Since 
the moment of inertia J^^i defined in Eq. (31), i.e., the 
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FIG. 4: The in-band B{E2) calculated by Eq. (46) for "=^Lu 
as a function of the rotational frequency. The mean-field pa- 
rameters are fixed and the same as in Fig. 3. Experimental 
data are taken from Ref. [16]. 



FIG. 5: The eigcnencrgy of the wobbling-like RPA solution 
as a function of the rotational frequency. The mean-field pa- 
rameters are fixed and the same as in Fig. 3. Experimental 
data are taken from Ref. [6]. 



so-called kinematic moment of inertia, contains the con- 
tribution of quasiproton alignment, ss i/'^rot (where i is 
the aligned angular momentum), the calculated J^^ de- 
creases as Wrot increases, if all the mean-field parameters 
are fixed. In contrast the measured inertia is almost con- 
stant against the change of the rotational frequency, sug- 
gesting at least one (or some) of the mean-field parame- 
ters is actually varying. The recently measured lifetimes 
of the TSD band [16] indicate that the experimental in- 
band B{E2), on the other hand, decreases rather rapidly, 
while the calculated in-band B{E2) is almost constant 
reflecting that the constant deformation parameters are 
used in this calculation. Therefore the B{E2) data also 
suggest that at least one (or some) of the deformation pa- 
rameters is changing as a function of Wrot ■ This behavior 
of the in-band B{E2) as well as out-of-band B{E2) will 
be discussed in more details in the following Sec. HIE. 



Since it is guaranteed in our formalism that the NG 
mode solution ti;„=NG = <^rot is precisely given by 
Eq. (18), we actually solve Eq. (35) to obtain the non-NG 
eigenenergy ujn and the amplitude (Xy^y (") : Xz^z (" ) ) by 
Eq. (34) , and then its forward and backward amplitudes 
are given by Eq. (23). As is mentioned in the previous 
section, it is important to use a sufficient model space 
in the RPA calculation for the consistent description of 
the NG mode (18). We have checked that the full use 
of the oscillator shells up to A^max = 12 is enough: Then 
the total number of proton and neutron two quasiparticle 
states (a, /3) is about 210,000 for the signature (— ) RPA 
calculation. Using the formulation explained in the pre- 
vious section, we have performed the Woods-Saxon RPA 
calculation, and found, just like in the previous Nilsson 
calculations [46-48], that the RPA modes that satisfy the 
requirements of the wobbling motion exist systematically 
in the TSD configurations. 



B. RPA calculations 

Once the nuclear mean-field is specified and the 
cranked quasiparticles, Eq. (11), are obtained, it is 
straightforward to solve the RPA equation, Eq. (25), or 
equivalently, Eq. (34) for the non-NG modes. It should 
be mentioned that the vacuum mean-field state of the 
RPA excitations is the lowest quasiproton excited state 
for the odd-Z nucleus ^^''Lu; such a state can be obtained 
by exchanging the HFB wave function and energy [27] of 
the excited quasiproton orbital to the opposite signature 
sector {U, V, E) {V, U, ~E). The details of the proce- 
dure to perform the RPA on such excited configurations 
are explained in Ref. [45], where the RPA excitations 
on the two quasineutron configuration (the Stockholm 
band) were treated. At the ground state (wrot = 0) any 
one quasiparticle configuration has the Kramers degener- 
acy so that it cannot be selected as a vacuum of the RPA, 
but at finite rotational frequency the degeneracy is lifted 
and then it is possible to perform the RPA calculation 
without any problems. 
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FIG. 6: The out-of-band to in-band B{E2) ratio as a function 
of the rotational frequency. The ratios associated with the 
7 — » 7 — 1 and 7 ^ 7 -f 1 out-of-band transitions are depicted 
as separate lines with the attached labels 7—1 and I + 1, 
respectively. The mean-field parameters are fixed and the 
same as in Fig. 3. Experimental data are taken from Ref. [16], 
in which only the 7^7—1 transitions are measured. 

As an example, the excitation energy of the wobbling- 
like RPA solution in ^^^Lu is compared with the ex- 
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FIG. 7: The calculated quadrupole amplitudes Q^^^Cn = wob) 
in Eq. (71) of the wobbling-like RPA solution as functions of 
the rotational frequency. 



perimcntal data in Fig. 5. In this calculation the con- 
stant mean-field parameters are used just like in Figs. 3 
and 4. The calculated excitation energy vanishes at a 
finite frequency, in this example about 0.5 MeV/ft, which 
is mainly due to the fact that many quasiparticle states 
come close to zero energy and the vacuum configuration 
changes as is seen in Fig. 2. The out-of-band B{E2) cal- 
culated by Eqs. (41) and (43) is an important quantity to 
identify the wobbling motion, which is as large as about 
100 Weisskopf units. In experiments the out-of-band to 
in-band B{E2) ratio, B{E2, 1 -> I - l)oui/B{E2, 1 
I — 2)in has been directly measured, and we compare the 
result of RPA calculation with experimental data for this 
ratio in Fig. 6. 

It should be noted that only the I ^ I — 1 transitions 
are measured for the out-of-band transitions, which sug- 
gests that the / ^ J-|- 1 out-of-band transitions are small 
and indicate the positive 7 shape. This characteristic be- 
havior of the rotor model also exist in the microscopic 
RPA calculation [42]. From Eqs. (41) and (43) the out- 
of-band B{E2) of the n-th RPA phonon excited state is 



B{E2;In^ I ±1y&c) 



i(Q^(n)TQf(n))'. (71) 



The electric quadrupole matrix elements Qy{n) and 
Qf{ri) have same sign and comparative magnitude for 
the wobbling-like solution (n = wob) for the positive 7 
shape, as they are depicted in Fig. 7, so that the 7 — > 7— 1 
transitions are much more enhanced than the 7^7 + 1 
transitions. In Fig. 6 the B{E2)'s of both transitions are 
shown; that of the 7 — > 7 -|- 1 transitions is smaller by an 
order of magnitude or more. 

The calculated energy of the wobbling excitation is 
100—200 keV smaller than the experimental data; this 
trend seems general and is similar in the previous Nils- 
son RPA calculations [46-48]. It is pointed out in 
Rcfs. [23, 24] that the particle-rotor coupling increases 
the excitation energy of the wobbling phonon consider- 
ably: Ziwp-wob ~ j/J^x in a simple limiting approxima- 
tion [23], where j is the angular momentum of the cou- 
pled quasiparticle so that /iwp-wob ~ 90 keV in the case 



of -'^^■^Lu. In our microscopic calculation the effect of 
the particle-rotor coupling is partly taken into account 
as a kind of blocking effect; the vacuum of the RPA ex- 
citation contains the odd quasiproton in the calculation. 
However, there exist explicit coupling effects between the 
RPA phonon and an odd quasiparticle. The underesti- 
mation of the wobbling excitation energy in odd-.^ Lu 
nuclei in the present and the previous RPA calculations 
may be due to this reason. More detailed theoretical 
investigations of the RPA-phonon-quasiparticle coupling 
calculation, see e.g. [77] , are necessary to draw a definite 
conclusion to this point. 

Although the calculated out-of-band to in-band B{E2) 
ratio in the present Woods-Saxon RPA calculation is 
about 10—20% larger than the previous Nilsson RPA cal- 
culation, it is still considerably smaller than the exper- 
imental data on average. Moreover, in contrast to the 
almost constant ratio of the experimental data, the cal- 
culated B{E2) ratio is decreasing as a function of the 
rotational frequency. Thus, the basic trend is similar to 
the previous Nilsson calculations, and was a most serious 
problem of the microscopic calculation compared with 
the macroscopic rotor model. It is, however, found that 
the problem of underestimation is due to the inconsistent 
definition of the triaxiality parameter "7" in the mean- 
field potential (either Nilsson or Woods-Saxon) and in the 
macroscopic model [55]. These characteristic features of 
the calculated and measured B{E2) will be discussed in 
more details in the following Sec. HIE. 
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FIG. 8: The three moments of inertia, and Jy%{n = wob) 
as functions of the rotational frequency. The mean-field pa- 
rameters are fixed and the same as in Fig. 3. 

Now let us discuss the results of calculation trans- 
formed into the PA frame. The three moments of inertia 
corresponding to the wobbling-like solution in Fig. 5 are 
shown in Fig. 8. The J'x inertia is the same as in Fig. 3. 
Other inertias. Jy^^ and J^°^ , are calculated according 
to Eq.(65) not to Eq.(38); the results of these two equa- 
tions are different for the general mean-field potential. 
We have found, however, that the actual difference be- 
tween the two definitions is small, within 1% at least in 
the present calculation; this suggests that the operators 
Fy z in Eq. (8) obtained with the present Woods-Saxon 
mean-field are not so much different from those propor- 
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tional to the quadmpole operators Qy^z in Eq- (14)- The 
calculated inertias J^^, are almost constant in this calcu- 
lation with fixed mean- field parameters, and Jy^ is much 
larger than . These basic features arc similar to the 
previous Nilsson RPA calculations [46]. We also show, 
in Fig. 9, the wobbling angles of the angular momentum 
vector in the PA frame (c.f., Eqs.(56), (59) and (60)), 
which are defined for the Ti-th RPA mode by [53, 54], 



tan^(n) = 



tan^(n) 



{JA 

^rl{n)+rl{n), 

{Jz)pA{n) ^ ry{n) 
{Jy)pA{n) TzinY 



(72) 



for which the small amplitude approximation requires 

tan0(n) « 9{n). From this figure the angle 9 is about 
25° (0.44 radian), and the small amplitude approxima- 
tion may be acceptable. Note that the quantities ry{n) 
and rz{n) describe the shape fluctuations, Eq. (59), and, 
at the same time, the angular momentum fluctuations, 
Eq. (60); ry{n) < rz{n) in this calculation and ry{n) de- 
creases rapidly as a function of Wrot (cf-, Fig. 7), and 
ry{n) <C rz{n) at higher frequency. This means that 
the angular momentum fluctuation of the y-direction is 
much larger than that of the 2-direction, and so the angle 
^ in Fig. 9. 
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FIG. 9: The wobbling angles of the angular momentum vector 
in the PA frame for the RPA wobbling excitation, 9{n = wob) 
and ■tp{n = wob), defined in Eq. (72), as functions of the 
rotational frequency. 

Finally, the calculated out-of-band B{M1) to in-band 
B{E2) ratio is compared with the experimental data in 
Fig. 10 as in the same way as the B{E2) ratio, but in 
the log scale. The B{M1) transition is calculated by 
Eqs. (42) and (44), i.e., 

B(Ml;/n^7±lvac) « i (/ij,(n) ± ^4"))', (73) 

where the RPA matrix elements of the Ml operator are 
given in the same way as those of E2 operator, 

IJ.y{n) = {[Xrr.,iHy]), IJ,z{n) = {[Xn,fiz])- (74) 
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FIG. 10: The out-of-band B{M1) to in-band B{E2) ratio as 
a function of the rotational frequency. Both tfic / ^ / — 1 and 
I —y I + 1 transitions are depicted as separate lines with the 
attached labels 7—1 and 7-1-1, respectively. Note that the 
ratios axe drawn in a log scale. The mean-field parameters 
are fixed and the same as in Fig. 3. Experimental data are 
taken from Rcf. [3], in which only the 7 — » 7 — 1 transitions 
are measured. 



The effective spin 51-factor gf^ = 0.6 gif®® is used [6, 22]; 
we have checked that the result with using gf^ = 0.7 g^J^^ 
does not change essentially. The 7^7—1 transitions 
are larger than the 7^7+1 transitions by factor two to 
four, in contrast to the case of B{E2) ratio. This is be- 
cause the lJ,y{n) amplitude is much larger than the Hzin) 
at all rotational frequencies in the present calculation. As 
shown in Fig. 10, the calculated ratio is about one order 
of magnitude larger than the experimental data, which 
is similar to that of Ref. [49]. The relative sign of the 
matrix elements between the E2 and A71 transitions is 
negative, which agrees with the experimentally measured 
^2/Ml mixing ratio [1, 6]. As for the B{M1) ratio, the 
measured values [3] are 0.0056 [/i^/e^fe^] at lower spins 
and 0.0098 at higher spins. The particle-rotor model cal- 
culations [1, 24] give about 0.015—0.02 at lower spins and 
about 0.01 at higher spins. In the previous Nilsson RPA 
calculation [46], the corresponding values are about 0.04 
and 0.015. Compared with these other calculations, our 
present value of B{M1) is too large, about 0.1, and al- 
most constant, while the other calculations predict that 
the B(M1) ratio decreases as a function of Wrot- We 
do not understand the reason why the present Woods- 
Saxon RPA calculation gives larger values especially at 
higher spin than the previous Nilsson RPA calculation. 
It should be mentioned that the absolute value of the 
B{M1) is small compared to the Weisskopf unit, but still 
it is much larger than the experimental data. This means 
that the RPA eigenmode operator corresponding to the 
wobbling excitation should be almost completely orthog- 
onal to the Ml operator (44), i.e., the various orbital 
contributions should cancel out, in order to obtain the 
agreement with the experimental data; the cancellation 
is not enough in the present calculation. Since our calcu- 
lation does not agree with experimental data for B{M1), 
we concentrate on the excitation energy and B(E2) in 
the following further investigations. 
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C. Dependence on mean-field parameters 

Until now we have shown the results of the micro- 
scopic calculations with fixed mean-field parameters at 
the reference values, {^2, J, Pi) = (0.42, 18°, 0.02) and 
= Ai, = 0.5 MeV. As it is explained in Sec. Ill A, we 
do not know very well about what arc the best values for 
these parameters. They might change depending on the 
rotational frequency, and then it must be clarified how 
their dependences are; for example, it is well-known that 
the proton and neutron pairing gaps should decrease as 
functions of the rotational frequency because of the Cori- 
olis anti-pairing effect. Instead of looking for what should 
be the values of these mean-field parameters, we here dis- 
cuss how the results of the present Woods-Saxon RPA 
calculation change when these mean-field parameters are 
varied in a reasonable range. 
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FIG. 11: Dependence of the eigenenergy (upper panel) or of 
the out-of-band to in-band B{E2) ratio (lower panel) on the 
pairing gap A = A-^ = Av for the wobbling-like RPA solution. 
The two curves correspond to the results at Wrot = 0.2 and 0.4 
MeV/ft. The mean- field parameters except the pairing gaps 
are the same as in Fig. 3. 

First, we show the dependence of the wobbling exci- 
tation energy and the B{E2) ratio on the pairing gap, 
A = Ajr = A^, i.e., the proton and neutron gaps being 
set equal for simplicity, in Fig. 11. In these figures, those 
at the rotational frequencies, Wrot = 0.2 and 0.4 MeV/H 
are shown. As it is clear the dependence of the wobbling 
excitation energy on the pairing gap is rather weak, and 
the very collective RPA solution exists even in the A = 



limit; at higher frequency cjrot = 0.4 MeV/fi. near the 
critical frequency of vanishing RPA energy, Wrot ~ 0.5 
MeV/h, the solution becomes unstable for larger pair- 
ing gaps. This is in contrast to the collective vibrational 
mode, and discussed [46, 48, 65] to be the character- 
istic property of the wobbling mode, which is of rota- 
tional origin. The dependence of the B(E2) ratio on 
A, which comes solely from the out-of-band B{E2), is 
slightly stronger, but again, smaller than that of the col- 
lective vibrational mode. The dependence of the crank- 
ing moment of inertia and the in-band B{E2) are well- 
known: The cranking moment of inertia decreases by 
about 27% (15%) at w^t = 0.2 (0.4) McV/ft as A in- 
creases from 0.1 to 0.8 McV for ^^'''Lu (not shown), while 
the in-band B{E2) does not depend on A essentially. 
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FIG. 12: Dependence of the eigenenergy (upper panel) or of 
the out-of-band to in-band B(E2) ratio (lower panel) on the 
deformation parameter [32 for the wobbling-likc RPA solution. 
The two curves correspond to the results at ojrot = 0.2 and 0.4 
MeV//i. The mean-field parameters except (32 are the same 
as in Fig. 3. 

Next, we show the dependence of the wobbling exci- 
tation energy and the B{E2) ratio on the magnitude of 
deformation (32 in Fig. 12. Other parameters are fixed 
at the reference values. The excitation energy slightly 
depends on P2 but is stable in the calculated range. The 
B{E2) ratio is almost independent of (32] this is quite 
reasonable in the macroscopic model, because both the 
in-band B{E2) and out-of-band B{E2) depend on the 
square of magnitude of deformation and it cancels when 
their ratio is taken. The in-band B[E2) increases mono- 
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tonically as (32 increases with keeping 7 (not shown), 
because the intrinsic quadrnpole moment about the x- 
axis increases. However, the cranking moment of inertia 
Jx is almost constant or shghtly decreases (not shown); 
this is because the core contribution increases but the 
quasiproton aHgnment contribution decreases and these 
two almost cancel. 
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FIG. 13: Dependence of the cigoncncrgy (upper panel) or of 
the out-of-band to in-band B{E2) ratio (lower panel) on the 
deformation parameter 7 for the wobbling-like RPA solution. 
The two curves correspond to the results at Wrot = 0.2 and 
0.4 MeV/R. The mean-field parameters except 7 are the same 
as in Fig. 3. 



Finally, we show the dependence of the wobbling exci- 
tation energy and the B{E2) ratio on the triaxiality pa- 
rameter 7 in Fig. 13. Again, the excitation energy slightly 
depends on 7 but is stable in the range, 15° < 7 < 30°. 
On the other hand, the B{E2) ratio increases rapidly 
as 7 increases. Both the in-band B{E2) and out-of- 
band B{E2) depend on the triaxiality in definite ways 
according to the macroscopic model. The result of the 
microscopic RPA calculation reproduces the characteris- 
tic features of the rotor model rather well, and this is 
precisely the reason that we can identify the obtained 
RPA solution as the wobbling excitation mode. The in- 
band B{E2) decreases monotonically as 7 increases (not 
shown, but see Sec. HIE). The cranking moment of iner- 
tia slightly increases (not shown); the core contribution 
decreases but the quasiproton alignment contribution in- 
creases, and the alignment contribution slightly wins in 
the present case. 



D. Calculation with varying pairing gaps 



The mean-field parameters are kept constants against 
the change of the rotational frequency in the results pre- 
sented above. However, at least the proton and neutron 
pairing gaps should decrease as functions of the rota- 
tional frequency due to the Coriolis anti-pairing effect. 
Therefore, it is desirable to determine and self- 
consistently. For this purpose, we fix the pairing force 
strength at the ground state with the pairing model space 
explained in Sec. Ill A, by 



(75) 



and use this strength for the TSD bands. Here Zi,-(co) 
is the experimental even-odd mass difference cvahiated 
by the third order difference formula [5], and {Pt) gi is 
the expectation value of the monopole pair transfer op- 
erator with respect to the normal deformed groimd state 
(wrot = 0), whose deformation is obtained by the Woods- 
Saxon Strutinsky calculation; i.e., /J2 = 0.217,7 = 
0,p4 = -0.002 for 163lu. 

The same Strutinsky calculation for the ground state, 
where the smooth pairing gap method is used with 
A = 13/VA MeV, gives the pairing gaps, ^7r,gr = 0.839 
MeV and Zi^,gr = 1.233 MeV for ^^^Ln. It 'should be 
noted, here, that the situation is puzzling: the even-odd 
mass differences for ^^'^Lu are Z\jr(co) = 1.243 MeV and 
A^{eo) = 0.868 MeV by using the 1995 mass data [78]. In 
the neighboring even-even nuclei, Ajrieo) = 1.180 MeV 
and Z\„(eo) = 1.167 MeV for ^^"^Yh, and Z\^(co) = 1.260 
MeV and Z\,.(eo) = 1.221 MeV for ^^-^Hf. Taking into 
account the blocking effect, the calculated pairing gaps 
{A^^gr = 0.839 MeV and Z\^,gr = 1.233 MeV) are reason- 
able in comparison with the even-odd mass differences of 
neighboring even-even nuclei. However, Z\7r(eo) in ^^^Lu 
is too large while Z\i,(eo) is too small compared with those 
in the neighboring even-even nuclei. This trend is not 
particular in ^''"^Lu, but it seems general in odd-proton 
nuclei in this rare earth region, for which we do not know 
the reason yet. In Ref. [60], we have used the neutron 
or proton pairing force strength which is an average of 
those for the neighboring even-even nuclei determined by 
the even-odd mass differences. In order to confirm that 
the results are not very different from those of Ref. [60] , 
we use, in the present work, the pairing force strengths 
determined by the even-odd mass differences of ^^^Lu. 

By using the pairing force strengths thus determined, 
we have done selfconsistent pairing calculation with 
the deformation parameters fixed at the reference value 
(/32,7,/34) = (0.42, 18°, 0.02). Then the proton and neu- 
tron pairing gaps vanish at Wrot ~ 0.42 and 0.43 MeV/fi, 
respectively, due to the alignments of quasiparticles. The 
abrupt pairing collapse is a drawback of the mean-field 
approximation, and we use the following phenomenolog- 
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FIG. 14: The proton and neutron pairing gaps used in the 
pairing varied calculation as functions of the rotational fre- 
quency for ^®^Lu. The deformation parameters used are 
(/32,7,/34) = (0.42, 18°, 0.02). 



ical parameterization at finite rotational frequency [79], 



^(Wrot) = ^0 X 



_^ _ 1 / Wrot Y 
2 V / 



1 / Uc 

2 \Wrot 



(76) 

where the parameters Aq = A{u)iot = 0) and lOc, such 
that Z\(aJrot = t^c) = 5^o> are determined by the selfcon- 
sistcnt pairing calculations. The resultant pairing gaps 
are shown in Fig. 14. 
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which is a great advantage over the previous Nilsson RPA 
model, where the cranking moment of inertia is overes- 
timated and the Strutinsky renormalization of the angu- 
lar momentum is necessary to reproduce the yrast TSD 
band. 



0.5 




EXP 



163 
71 



Lu 



0.3 0.4 
h«,ot [MeV] 



0.5 



0.6 



FIG. 16: The cigcncncrgy of the wobbling-liko solution as a 
function of the rotational frequency. The deformation param- 
eters are fixed but the pairing gaps in Fig. 14 are used. 
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FIG. 17: The three moments of inertia, J^x and Jy%{n = wob) 
as functions of the rotational frequency. The deformation 
parameters are fixed but the pairing gaps in Fig. 14 are used. 
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FIG. 15: The cranking moment of inertia as a function of the 
rotational frequency. The deformation parameters are fixed 
but the pairing gaps in Fig. 14 are used. 

We have done the semi-phenomenological pairing self- 
consistent calculation by using the pairing gaps in Fig. 14 
and the constant deformation parameters at the reference 
values, (/32,7,/34) = (0.42, 18°, 0.02). The result for the 
cranking moment of inertia is shown in Fig. 15. Com- 
paring with the constant pairing calculation in Fig. 3, 
the agreement between the experimental data and the 
calculation is much better. This means that the present 
Woods-Saxon RPA model can describe both the yrast 
TSD band and the wobbling excitation at the same time, 



The result of the excitation energy of the RPA wob- 
bling mode is compared with the experimental data in 
Fig. 16. The behavior of the excitation energy is not 
very different from the constant pairing calculation in 
Fig. 5, but the overall energy is slightly lower and about 
200 keV smaller than the experimental data. The three 
moments of inertia in the PA frame calculated with vary- 
ing pairing gaps are shown in Fig. 17. In contrast to the 
corresponding constant pairing calculation in Fig. 8, the 
inertias Jy^^ acquire sizable dependence on the rotational 
frequency, but the basic features such as Jy^ ^ J^^ 
are almost the same. As for the B{E2)^s, the in-band 
B{E2) and the out-of-band to in-band B{E2) ratio are 
shown as solid lines in Figs. 18 and 19. As it is dis- 
cussed in Sec. Ill C, the in-band B{E2) does not depend 
on the pairing gaps, so the result is almost the same as 
that with the fixed mean-field parameters in Fig. 4. The 
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B{E2) ratio is slightly different from that in Fig. 6, but 
the difference is small. Thus, the calculation with varying 
pairing gaps is not helpful to fill the gap for the discrep- 
ancy between the calculation and the experimental data 
of B{E2ys, which is considered in the next subsection. 
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FIG. 18: The in-band B{E2) defined for ^'^'''Lu as a function of 
the rotational frequency (solid line). The deformation param- 
eters are fixed but the pairing gaps in Fig. 14 are used. The 
dashed line is the result of calculation with further changing 
the 7 parameter Unearly from 7 = 18° at Wrot = 0.2 MeV/h 
to 7 = 30° at Wrot = 0.4 MeV/h. 



interband E2 transitions. It is important not only to 
identify the wobbling excitation but also to deduce the 
triaxiality parameter 7; according to the rotor model, the 
out-of-band to in-band B{E2) ratio directly reflects the 
triaxiality of the quadrupole deformation. Here it should 
be noted that there are various different definitions of the; 
triaxiality parameter 7, and one should be careful which 
definition is used for quantitative discussions [55]. 

The parameter 7 in the moan- field potential, Eq. (5), is 
one of such definitions, called 7(pot:WS) in Ref. [55], but 
it is different from the one specified by the quadrupole 
moments, which is called 7(den) and defined by 



(77) 



where Q is the total quadrupole moment. In terms of 
this 7 = 7(den), the quantities Uy^zR^ in Eq. (50) can be 

expressed as UyE? = —yj Qsm{'y + 60°) and azR"^ = 

~\j W '3sin7, and then by using Eq. (67) with (59) the 

out-of-band B{E2) in Eq. (71) is represented as in the 
same form as the macroscopic rotor expression [5], 
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FIG. 19: The out-of-band to in-band B{E2) ratio for "^Lu as 
a function of the rotational frequency (solid line). The defor- 
mation parameters are fixed but the pairing gaps in Fig. 14 
arc used. The dashed line is the result of calculation with 
further changing the 7 parameter linearly from 7 = 18° at 
cjrot = 0.2 McV/ft to 7 = 30° at Wrot = 0.4 MeV/fi. The 
B{E2) ratio with only I —^ I — 1 transitions is shown. 
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Here {J^) « 7, Q^^ = (eZ/^)Qj;.^, and c; 



Wy{n) = 1/jfin) - 1/J,, Wzin) = l/jf{n) - 

(79) 
1 as well 

as (T = + in Eq.(68) arc further assumed. As it is dis- 
cussed in Sec. II C, the condition = 1 and a = + 
can be regarded as a requirement for the wobbling-like 
RPA solution: In the present and the previous calcula- 
tions this condition is satisfied if the used model space 
is large enough to guarantee the exact NG mode decou- 
pling; the deviation is within 2% in the present calcula- 
tions. Using the same 7 = 7 (den) and the assumption 
= {eZ/A)Q2u, the in-band B{E2) in Eq.(46) can be 
also expressed as 



E. B{^E2) and triaxial deformation 

The out-of-band B{E2) from the excited band to the 
yrast TSD band is the crucial quantity to identify the nu- 
clear wobbling motion. The recent lifetime measurement 
of the excited and yrast TSD bands [15, 16] revealed that 
the out-of-band B{E2) is as large as about 100 Weisskopf 
units, which is much larger than the typical coUectivc- 
vibrational E2 transitions, and is possibly the largest 



BiE2-J^I±2)^±-(e'-Q 



cos2(7-F30°). (80) 



In this way, the out-of-band to in-band B{E2) ratio does 
not depend on the total quadrupole moment Q (and so 
not on /32 and (34) in a good approximation, and the de- 
creasing trend of the calculated B{E2) ratio as a function 
of the rotational frequency can be naturally understood 
by the 1/7 factor in Eq. (78). 

It should be emphasized that the 7 parameter that ap- 
pears in Eqs. (78) and (80), which are essentially the 
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formula employed in the particle-rotor model calcula- 
tions [22-25] (note that the effect of the odd particle 
is negligible for B{E2)'s), is 7 = 7(den) in Eq. (77) 
and not 7(pot:WS) in Eq. (5). As it is demonstrated 
in Ref. [55], they are rather different at a given shape 
especially for larger quadrupole deformations; for exam- 
ple, 7(den) = 20° corresponds to 7(pot:WS) » 30°, for 
the TSD shape of the Lu and Hf region (sec Fig. 4 of 
Ref. [55]). It was discussed in the particle- rotor model 
calculations [22, 23] that 7(den) ft! 20° are necessary in 
order to reproduce the experimental B{E2) ratio on av- 
erage. This means that one has to use 7 » 30° in the 
Woods-Saxon potential in Eq.(5). Since the calculated 
B{E2) ratio increases as a function of the parameter 7 
as is discussed in Sec. Ill C, we obtain similar amount of 
agreement to that in the particle-rotor model if the in- 
creased value 7 ft 30° is used: This is demonstrated for 
the Nilsson RPA calculation in Ref. [55]. 

Although the average magnitude is improved by in- 
creasing the triaxiality of the mean-field, there is still con- 
siderable disagreement between the calculated and mea- 
sured B{E2) ratios: Their rotational frequency depen- 
dences arc different. This suggest that some parameter 
of the mean-field should change as a function of the ro- 
tational frequency. Here it should be mentioned that the 
measured in-band B{E2) also depend on the rotational 
frequency; it decreases as Wrot increases. This trend can 
be understood either as a result of decreasing or 
as a result of increasing 7 (or both of them). However, 
the B{E2) ratio is almost independent on fi2 [Pa)- There- 
fore, a simplest way to understand both B{E2) data is 
to increase 7 as a function of Wrot- In the present work, 
we try to reproduce the B{E2) data by varying the 7 
parameter linearly from 7 = 18°, i.e., 7(den) ft 12°, at 
Wrot = 0.2 MeV/fi to 7 = 30°, i.e., 7(den) ft 22°, at 
^Tot = 0.4 MeV/h. The results are shown as dashed 
lines in Figs. 18 and 19, where rather good agreements 
for both the in-band B{E2) and the B{E2) ratio are ob- 
tained. Considering that the various Nilsson Strutinsky 
calculations [14], including ours, predict that the triax- 
iality parameter does not change so largely, we need to 
look for more reliable methods to determine the mean- 
field parameters if this amount of change of the triaxial 
deformation is really true. 



F. Systematic calculations in Lu and Hf nuclei 

Until now we have considered only the nucleus ^^^Lu as 
a best example. However, the wobbling excitations have 
been identified in several neighboring isotopes, -"^^-^Lu [10], 
^^^Lu [8], and ^^''Lu [9]. Moreover, as it was pre- 
dicted [11, 12, 14], the many possible candidates of the 
TSD bands have been observed in even-even Hf iso- 
topes [17-21]. Therefore, it is interesting to know that 
the existence of the wobbling excitation is specific for 
the observed Lu isotopes or a general phenomenon. In 
our microscopic RPA formalism, the wobbling mode is a 



collective excitation mode like the low-lying collective vi- 
brations, the (3- or 7-vibrations, which are well-known to 
exist systematically near the ground states in the wide 
range of the periodic table. It can be inspected that 
the wobbling-like RPA solutions will be found gener- 
ally in TSD nuclei. Thus, we have performed system- 
atic RPA calculations for several Lu and Hf isotopes 
with A?^ = 90 — 96. Although the mean-field parameters 
must be different in each isotope, we cannot determine 
them so that we have used for all nuclei the fixed ref- 
erence values that are used for ^^^Lu in Sec. IIIB, i.e., 
(/32,7,/34) = (0.42, 18°, 0.02) and = = 0.5 MeV. 
The results are shown in the following four figures: The 
wobbling excitation energies of the Lu and Hf isotopes 
are shown in Figs. 20 and 22, and the out-of-band to in- 
band B{E2) ratios of the Lu and Hf isotopes in Figs. 21 
and 23, respectively. 
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FIG. 20: The cigcnoncrgics of the wobbling-like RPA solu- 
tions in Lu isotopes (Z = 71) with =90, 92, 94, and 96 as 
functions of the rotational frequency. The mean-field param- 
eters are fixed and the same as in Fig. 3. 
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FIG. 21: The out-of-band to in-band B{E2) ratios in Lu iso- 
topes with A'' =90, 92, 94, and 96 as functions of the rotational 
frequency. The mean-field parameters are fixed and the same 
as in Fig. 3. 

As it is clear in these figures, the obtained wobbling ex- 
citation energies are rather similar for the isotopes, and 

are about 180 (220) keV in the Lu (Hf) isotopes, with 
the exception of A^ = 96 isotope, whose excitation energy 
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FIG. 22: The eigenenergies of the wobbhng-hke RPA solu- 
tions in Hf isotopes {Z = 72) with N =90, 92, 94, and 96 as 
functions of the rotational frequency. The mean-field param- 
eters are fixed and the same as in Fig. 3. 



occupied. On the other hand the vacuum states are zero 
quasiproton states in the Hf isotopes, and the lowest two 
quasiprotons align at Lo^ot ~ 0.23 MeV/h with relatively 
large interactions. As it is discussed in Refs. [46, 48] , the 
proton quasiparticle alignment is an important factor to 
stabilize the RPA wobbling excitation for the positive 7 
shape, so that the wobbling-like solutions appear only 
after the two quasiproton alignment in the Hf isotopes. 

It must be emphasized that the collective RPA excita- 
tion exists irrespective of the proton number being odd 
or even in our microscopic RPA formalism. Therefore, 
the wobbling excitations should be identified also in the 
even-even Hf nuclei. It is mysterious for us that there has 
been no sign of them in any Hf isotopes: This may sug- 
gest that some important elements may be still missing, 
which we do not know yet. 
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FIG. 23: The out-of-band to in-band B{E2) ratios in Hf iso- 
topes with N =90, 92, 94, and 96 as functions of the rotational 
frequency. The mean-field parameters are fixed and the same 
as in Fig. 3. 



increases by about 100—200 keV in Wrot > 0.3 MeV/h. 
The reason why the excitation energy increases in the 
= 96 isotope is that an additional two quasineutron 
alignment occurs gradually around that rotational fre- 
quency in this particularly chosen mean-field parameters. 
We do not think that it is very meaningful; in fact the 
observed wobbling excitation energy in ^^^Lu is very sim- 
ilar to that in ^^^Lu, and we have checked that the slight 
modification of the deformation parameters changes this 
trend of increasing energy of the N = 96 isotopes. The 
small kinks of = 90 isotopes at Wrot ~ 0.32 — 0.35 
MeV/h are caused by a sharp crossing of two quasineu- 
trons, whose orbital has rather small alignment and af- 
fects the vacuum very little. This is also a kind of artifact 
due to the mean-field parameters used presently, since 
the sharp crossing is not observed in experiment for the 
-'^^-'^Lu isotope. The B{E2) ratios in the Lu and Hf iso- 
topes are also very similar, again except for A^ = 96 iso- 
topes, whose B{E2) reduce slightly at higher rotational 
frequency due to the quasineutron alignment. In the Lu 
isotopes, the odd high-j [Nqsc = 6) quasiproton is always 



IV. SUMMARY 

We have investigated the nuclear wobbling motion in 
the Lu and Hf region by using the microscopic frame- 
work but in the small amplitude approximation, i.e., the 
cranked mean-field and the random phase approximation 
(RPA). The concept of the wobbling motion is intimately 
connected to the rotor model: The body-fixed frame or 
the principal-axis (PA) frame, where the degrees of free- 
dom corresponding to the collective rotations are elimi- 
nated, plays a key role. In the conventional dcsc;ription, 
i.e., in the uniformly rotating (UR) frame of the crank- 
ing prescription, in which no concept of the body-fixed 
frame is introduced, the wobbling excitation is just like 
an usual collective vibration; only after transforming into 
the PA frame it can be interpreted as a motion of the 
angular momentum vector wobbling about the main ro- 
tation axis. Therefore, we reviewed the microscopic RPA 
framework [32], where the PA frame is introduced, from 
a slightly general view point [42]: The original frame- 
work by Marshalek [32] assumes that the deformation of 
the mean-field is purely of quadrupole type, but such a 
limitation should be lifted for general mean-field poten- 
tials. By introducing the symmetry-preserving residual 
interaction associated with a given mean-field potential, 
it is shown that the original Marshalek formulation can 
be recovered with a slight modification. Furthermore, 
the strength of this residual interaction is uniquely de- 
termined for a given mean-field so that there is no new 
parameter in the calculation of the RPA step, i.e., the 
theory always has the consistency between the mean-field 
and the residual interaction. 

As a mean-field, in the present work, we have used the 
new Woods-Saxon potential that gives correct density 
distributions for both protons and neutrons. It is be- 
lieved to give better descriptions than, e.g., the Nilsson 
potential, which was employed in the previous works [46- 
48]. It has been found that the RPA solutions, which can 
be interpreted as wobbling excitations, exist systemati- 
cally in the Lu and Hf region, if the proper mean-field 
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parameters arc adopted, i.e., the pronounced triaxial de- 
formation or the triaxially superdcformed (TSD) shape. 
The general features of the calculated RPA solutions are 
rather similar to those obtained in the previous works em- 
ploying the Nilsson potential, except for the B{M1) ratio 
that is too large, and the stability against the change 
of the mean-field parameters are demonstrated in de- 
tails. It should be emphasized that the obtained RPA 
solutions satisfy the properties expected from the rotor 
model; especially the large out-of-band B{E2) is repro- 
duced if the same triaxial deformation as the particle- 
rotor model analysis [22-25] is used. All these results 
confirm the wobbling excitations in the TSD bands from 
the microscopic view point. 

We must admit that the results of our calculations do 
not agree very well with the experimental data of the 
Lu isotopes. The excitation energies of the wobbling ex- 
citation are underestimated about 150—200 keV, which 
may suggest that the particle-rotor coupling effect is not 
fully taken into account. We obtain the similar excitation 
energies for wobbling mode in the even-even Hf nuclei. 
Therefore, it is crucial to observe the wobbling excita- 
tions in even-even nuclei to draw a definite conclusion 
for whether the effect of the particle-rotor coupling is es- 



sential for the wobbling excitation observed in the odd-Z 
Lu isotopes. As for the out-of-band to in-band B{E2) 
ratio, the calculated results decrease as functions of the 
rotational frequency just like in the rotor model, if the 
deformation parameters are kept constants. The mea- 
sured B{E2) ratio in ^^'"'Lu is almost constant, and one 
has to increase the triaxiality parameter considerably as 
a function of the rotational frequency in order to repro- 
duce this behavior. However, the Nilsson Strutinsky cal- 
culation gives almost constant triaxial deformation and 
the obtained triaxiality is too small to account for the 
average magnitude of the B(E2) ratio. The better mi- 
croscopic description of the wobbling motion is certainly 
a challenge to the existing microscopic theory. 
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